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Structural Disorder Induced Polaron Formation and Magnetic Scattering in the 
Disordered Holstein-Double Exchange Model 

Pinaki Majumdar and Sanjeev Kumar 

Harish- Chandra Research Institute, 
Chhatnag Road, Jhusi, Allahabad 211 019, India 

In this paper we present results on the disordered Holstein-Double Exchange model, explicitly 
in three dimension and 'metallic' densities, obtained by using a recently developed Monte Carlo 
approach. Following up on our earlier paper, cond-mat 0406085, here we provide a detailed 
microscopic picture of the thermally driven metal-insulator transition (MIT) that arises close 
to the ferromagnet to paramagnet transition in this problem. This paper is focused mainly 
on the 'diagnostics', clarifying the origin of the effective disorder that drives the MIT in this 
system. To that effect, we provide results on the thermal evolution of the distributions of 
(i) lattice distortions, (ii) the net 'structural disorder' and (iii) the 'hopping disorder' arising 
from spin randomness feeding back through the Hunds coupling. We suggest a phenomenology 
for the thermally driven MIT, viewing it as an 'Anderson-Holstein' transition. 

KEYWORDS: colossal magnetoresistance manganites, double exchange, lattice polarons, electron locali- 
sation 
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1. Introduction 

The Holstein-Double Exchange (H-DE) model has 
been the focus of recent attention because it can help 
clarify the combined effect of strong electron-phonon 
(EP) interaction and strong electron-spin coupling that 
operates in the manganites, Lai-^Ca^MnOa, say. 1 The 
detailed microscopic model for manganites is rather com- 
plicated. Let us start by describing this model, and argue 
how successive reduction to the disordered H-DE model 
retains some of the essential physics. 

The comprehensive tight binding model for perovskite 
manganites, neglecting oxygen orbitals, etc, is 

a/3 

H = Yl ^fcLaCjp* + e i n i - J H Y S i-°i 
(ij)a i i 

— A ^ Qi-fj + -ffstiff + i?ph-dyn + -ffHubb (1) 

i 

The tf.® are hopping between doubly degenerate Mn e g 
levels at neighbouring sites, refers to weak substitu- 
tional disorder, and Jr is the strong Hunds coupling be- 
tween the e g electrons and the 5 = 3/2 'core spin' fag 
electrons). The 'orbital moment', r,-, of the e g electrons is 
Jahn- Teller coupled to the octahedral distortion param- 
eter Qi, through the EP coupling A. -ff st iff is the stiffness 
of the lattice. In its simplest form it is ~ (if/2) XL Qf; 
but in reality is of 'cooperative' character, i.e, involves 
phonon degrees of freedom in more than one octahedra. 
The phonons are quantum variables, with intrinsic dy- 
namics arising from ff p h-dyn ~ P 2 /(2Mq), where P is 
the momentum and Mq the mass of the relevant oscil- 
lator. Finally, the x = state in the manganites is a 
Mott insulator, arising (partly) from large on site 'inter- 
orbital' Hubbard repulsion, i?Hubb ~ U'J^i a ^f3 n ia n i/3- 
The electron-phonon, electron-spin and Hubbard in- 
teractions are all large 3> t, where t is the typical hopping 
scale in the problem, and therefore beyond the range of 



perturbation theory. The detailed model, unfortunately, 
is far too complex for the present methods of many body 
theory, if it were to be handled in a realistic three di- 
mensional situation. However, it is possible to simplify 
the model somewhat, recognising that (i) at large dop- 
ing of the Mott insulator, e.g, x ~ 0.3 — 0.4, the Hubbard 
interaction probably does not have a qualitative effect, 
(ii) the phonons are in the adiabatic regime, with typical 
frequency w p h <C t and, as a first approximation, we can 
explore the adiabatic limit w p h = 0, and (iii) the S = 3/2 
may be approximated as a 'classical' spin. 

This still leaves us with the problem of (disordered) 
electrons strongly coupled to (classical) phonon and spin 
degrees of freedom. The qualitative effects in such a sys- 
tem include (i) magnetic order, typically ferromagnetism, 
arising from Hunds coupling and electron delocalisation, 
(ii) the possibility of polaron formation if the EP cou- 
pling is sufficiently large compared to the kinetic scale, 
and (iii) the phase competetion between different kinds 
of long range order. These complexities, observed exper- 
imentally in the manganites, can be recovered already at 
the level of the Jahn- Teller-Double Exchange (JT-DE) 
model, neglecting 7J ph _ dyn and izHubb- 

The JT-DE model was explored by Millis et al., 2 
within dynamical mean field theory (DMFT), and later 
by Dagotto and coworkers 3 using real space Monte Carlo 
(MC) tecniques. The DMFT study could demonstrate a 
thermally driven metal- insulator transition (MIT), as ob- 
served in Lan.7Cao.3Mn03, although only at x = 0.5. The 
MC studies, on the other hand, have gone a long way in 
clarifying the various phases in the JT-DE model, but 
the severe finite size handicap has prevented estimate of 
transport properties. 

If we tentatively accept that the essential physics in the 
manganites arises from the general interplay of EP cou- 
pling and double exchange, then the JT-DE model can 
be further reduced to a one band version: the Holstein- 
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Double Exchange model, defined further on. Indeed, it 
was studied early on, 4 simultaneously with the JT-DE 
work of Millis et al, within a mean field approximation, 
to clarify the effect of EP coupling on the ferromagnetic 
T c . That work did not address transport properties. More 
recently the H-DE model has been studied via many 
body CPA approximations 5 and through direct Monte 
Carlo simulation. 6 The MC study does find a thermally 
driven MIT but only at low carrier density. In addition, 
all these studies neglect the effect of quenched disorder 
which should have a dramatic impact 7 in a system with 
strong EP coupling. 

The effect of disorder in 'manganite models' have in 
fact been discussed, 8 ' 9 but the focus has been on its 
impact on phase competetion and bicriticality. Disorder 
can lead to nanoscale cluster coexistence near a first or- 
der phase boundary 10 and much of the disorder related 
theory has concentrated on modelling this phenomena. 
While this is a vital issue, there is a completely indepen- 
dent effect of disorder that shows up in systems with 
strong EP coupling. Disorder can induce polaron for- 
mation, even below the polaronic threshold for 'clean' 
systems, and drastically modify the residual resistivity, 
the optical spectral weight, and in fact the entire trans- 
port response. Manganite data testifying to these effects 
abound in the literature, but there has been no theoret- 
ical effort to understand these phenomena. 

2. Experiments in the Manganites 

In the manganites, AjrA'i-^MnOs say, the physics is 
controlled by varying x, or the mean A-A' site cation 
radius ta, or the cation disorder a a, arising from size 
mismatch. Since the parameters and a a themselves 
depend on x, it is best to focus on data at a fixed dop- 
ing level. In this spirit, experimenters have used chemi- 
cal variation at fixed x, varying both rA (which controls 
the electronic bandwidth) and a a- 11,12 Pressure experi- 
ments, on the other hand involve variations in the band- 
width, keeping the disorder constant. 13, 14 Finally, the 
effect of disorder has been explored, keeping x and ta 
fixed, varying only a a- 15,16 

The main conclusion of these studies is that (i) large 
ta systems, Lai-^Sr^MnOa, are 'canonical DE magnets', 
with a saturated ferromagnetic ground state, 'metal- 
lic' resistivity at all temperature, and modest mag- 
netoresistance (MR) near T c , (ii) reducing ta, as in 
Lai_ x Ca x Mn03, maintains a saturated ferromagnetic 
ground state, but there is a thermally driven MIT near 
T c , with associated colossal magnctoresistance, and a 
pseudogap in density of states (DOS), (iii) further re- 
duction in rA as in Lai-x-yPryCa^MnOs, leads to a 
non saturated ferromagnetic ground state, mixed phase 
tendency, nanoscale clusters (depending on disorder), a 
field driven insulator-metal transition at zero tempera- 
ture, and large MR over a wide temperature range. All 
these occur with only a few percent change in rA ■ 

These data define the problem for a 'global' descrip- 
tion of the effects seen in the manganites. With small 
variation in electronic bandwidth (or inversely the EP 
coupling), ~ 10%, a theory should be able to reproduce 
the suppression in T c , the dramatic changes in transport 
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character, reduction of low frequency optical spectral 
weight, sharply enhanced residual resistivity, and emer- 
gence of a pseudogap features in the DOS. Simultane- 
ously it should reproduce similar effects if the quenched 
disorder were increased, remaining at fixed EP coupling 
and bandwidth. 

In our earlier paper 17 we had solved the disordered 
Holstein-Double Exchange (d-H-DE) model and qualita- 
tively reproduced all the trends above (except magnetic 
phase competetion, which requires additional antifcrro- 
magnetic coupling). In this paper we focus on a generic 
point in parameter space, where there is a thermally 
driven MIT, and try to provide a detailed microscopic 
picture for the transition. 

3. Model and Method 

The d-H-DE model, with classical spins and phonons, 
is defined by: 

h = -t 4a c j<? + J2( e * ~ ~ Jh St - at 

- A^Tn^ + ^-^Tx 2 (2) 

i i 

Here the t are nearest neighbour hopping on a sim- 
ple cubic lattice, ej is the quenched binary disorder, with 
value ±A, and Jh is the Hunds coupling. A is the EP 
interaction, coupling electron density to the local distor- 
tion Xi, and K is the stiffness of the lattice. The basic 
parameters in the problem are A/i, A/i, electron density 
n, and temperature T. We work in the limit Jujt — > oo. 
We also set K = 1, and measure energy, frequency, T, 
etc, in units of t. 

For Jn/t — > oo, the Hunds coupling acts as a con- 
straint, orienting the electron spin at a site parallel to 
the core spin (and projecting out the anti-parallel com- 
ponent) leading to the following spinless fermion model: 

H = + ^c) + $>^ + y]Tx 2 

(ij) i » 

Uj = i(cos|cos^+i sin I sin ^e^"^)) 
J y 2 2 2 2 

& = £j - IJ- ~ (3) 

The 7, 7^ are spinless fermion operators, corresponding 
to the 'parallel' spin projection of the original electrons. 
£i is the net 'potential' seen by the electrons, while the 
hopping 'disorder' is controlled by 6i and fa, the polar 
and azimuthal angles, respectively, of S^. 

The problem involves strong coupling and disorder, 
and needs to handle thermal fluctuations. We use a re- 
cently developed Monte Carlo technique 18 based on a 
"travelling cluster approximation" (TCA) to anneal the 
classical variables. This allows us to use system sizes 
N ~ 10 3 and address the MIT in the model. 

After equilibriating the phonon and spin degrees of 
freedom for a specified set of parameters, disorder reali- 
sation, and T, we compute the following: (i) transport 
properties, adapting a scheme 19 which now uses elec- 
tronic eigenfunctions in the annealed background, and 
Fermi factors, for T ^ 0, (ii) DOS, (iii) the thermally 
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Fig. 1. Panel (a), shows the variation in the magnetisation, m(T), 
the mean magnitude of the hopping amplitude, t av (T), and the 
r.m.s fluctuation in the hopping, 8t(T). Panel (6). shows the 
scaled resistivity, p(T), the scaled variance of the lattice dis- 
tortion, Sx 2 , and the variance of the effective disorder Srj 1 . 
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Fig. 3. The temperature dependence of the 'd.c' resistivity at 
three system sizes, L = 8, 10, 12, with N = L 3 . 
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Fig. 2. (a) The low energy density of states, and (a) the low 
frequency optical conductivity, with varying temperature. 



averaged lattice distortion, as well as (iv) the distribu- 
tion of net 'structural disorder' and hopping disorder, 
and (v) spatially resolved information on density dis- 
tribution, n r , and the magnetic correlation, / ma g(r) = 
^ r ,(S r .S r '}, where r, r' are nearest neighbours. 

4. Results 

In our earlier, main paper, 17 we have provided the 
n — T phase diagram in the clean limit for various A, 
while the transport, spectral and optical properties were 
shown for n = 0.3 varying the EP coupling and disor- 
der. We demonstrated a thermally driven MIT over a 
wide parameter regime in A — A but did not have room 
to provide a detailed microscopic picture of the 'ther- 
mal disorder' that drives the MIT. This paper focuses 
on a single, generic, point in parameter space n = 0.3, 
A = 2.0, A = 0.6 and tracks the T dependence of the 
distribution of effective disorder as well as the spatially 
inhomogeneous underlying state. 

If the thermally driven MIT in the manganites is a 
"localisation transition" , rather than simple "band split- 
ting", a real space understanding of the phenomena 
should be of crucial importance. 

Figure 1 (b) shows the resistivity at our chosen param- 
eter point. Figure 1(a) shows the variation in magnetisa- 
tion, m(T), the spin disorder induced suppression of the 
average, t av — (\Uj\), and the r.m.s fluctuation St(T) in 
the hopping distribution. 
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Fig. 4. Size dependence in the low frequency density of states, 
showing the stability of the pseudogap feature. System sizes, N = 
L 3 with L = 8, 10, 12, as in Fig.3. 



The evolution of the mean square lattice distortion, 
as well as the effective structural disorder arising from 
ej — Xxi are shown in panel (b). The full distribution of 
these disorder is shown in Fig. 5. The kink like feature in 
p(T) is a consequence of a discrete set of points sampled 
in T. Denser sampling reveals a more continuous change 
in character. 

Apart from the MIT, whose detailed description is our 
primary focus in this paper, Fig. 2(a) shows the thermal 
evolution of the low energy DOS, while Fig. 2(b) shows 
the low frequency optical conductivity. The two features 
to note are (i) the non monotonic T dependence of the 
low energy DOS (the 'dip' deepens initially with increas- 
ing T and then fills up) and (ii) the strongly non Drude 
nature of the optical response even at T = 0. The resis- 
tivity and DOS had been shown in our earlier paper. 17 In 
what follows we want to focus on the (annealed) disorder 
that is responsible for these physical effects. The finite 
size effects in this problem are rather weak, as borne out 
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Fig. 5. Panel (a): P(x) the distribution of lattice distortions, av- 
eraged thermally and over the system, (a) the distribution of 
the effective potential, r\i = e; — Arr;, seen by the electrons, and 
(c) the distribution in the magnitude of the hopping, essentially 
the nearest neighbour correlation sj\ + Si.Sj. 



by the size dependence of the resistivity in Fig. 3, and 
the density of states, Fig. 4. 

Effects for T — > 0: Figure 1(b) shows that there are lat- 
tice distortions of fairly large magnitude even at T = 0, in 
the fully polarised ferromagnetic state. An understand- 
ing of this can be obtained by studying the disordered 
spinlcss Holstein model 20 where we observe that close 
to a (collective) polaronic instability, even weak disorder 
can induce large lattice distortions and localise a fraction 
of charge carriers. This phase differs from a polaronic in- 
sulator in that there are still extended states that sur- 
vive close to the Fermi level. The DOS shows signature 
of this partial polaronic trapping (i.e a pseudogap), and 
the optical response has the non Drude character typical 
of electrons in a strongly disordered background. 

Figure 5 shows the actual distributions, of lattice dis- 
tortions, Xi, the net potential rji — — Xx i} and the 
magnitude of hopping, \Uj\. The quenched disorder is bi- 
nary with value ±0.6. Starting with panel (c), the hop- 
ping distribution P(\t\) tends to a 5 function as T — > 
since all spins are aligned. In this regime the physics 
is controlled purely by the effective structural disorder. 
The P(x) distribution is bimodal, with significant weight 
near x ~ X/K = 2, which would arise for strongly lo- 
calised particles. There is also a peak related to the orig- 
inal clean Fermi liquid, surviving near the origin. 

We can crudely estimate the number of strongly lo- 
calised carriers in terms of dxP(x) where x m ; n = 1.2 
is the (approximate) upper limit of distribution at weak 
disorder, A = 0.2, Fig. 7. This area is roughly 0.2 - 0.25 
at A = 0.6. Since the electron density itself is n = 0.3 this 
suggests that a large fraction of carriers have collapsed 
into "polaronic" states. The density profile at T = 0.01, 
Fig. 6, top right, confirms this picture. That not all par- 
ticles are localised is borne out by the residual resistivity, 
p(0). We have checked that large L extrapolation of our 
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Fig. 6. The view on the top surface of a 10 system. The row 
above shows the thermally averaged density profile, n r , at 
T = 0.01,0.08,0.30 (right to left). The row below shows the 
nearest neighbour averaged magnetic correlation: fmag(r) = 
S r /(S r .S r /), for the same sequence of temperatures. 
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Fig. 7. The distribution P(x, A) at T = 0.01 showing the evolu- 
tion of the polaronic peak with increasing disorder. Inset: vari- 
ation in the density of localised carriers, n po ;, with increasing 
disorder, staying at A = 2.0 and fixed electron density, n = 0.3 
We have used n po ;(A) = dxP(x,A) with x m ; n = 1.2. 



results (verified using N = 6 3 — 12 3 ) still leads to finite 
P(0). 

The P(rj) distribution, at T = 0, is more complex than 
P(x), and is more relevant for understanding the 'land- 
scape' in which the electrons move. As a starting point we 
can imagine three kinds of sites: (i) sites with ej = ±A, 
with weak distortions and no 'overdensity', these sites 
would be continuation of the clean Fermi liquid, call them 
vfl, (ii) sites with m = — A and some large distortion 
(call these r\~) and (iii) sites with ej = +A and moder- 
ately large distortions (?7 + ). 

The xi for a 'site localised' electron is ~ X/K = 2, 
so the n~ sites would have magnitude w —A — X 2 / K = 
—4.6. In fact the 'polarons' are not quite site localised, so 
r\~ would be somewhat less: consistent with the left peak 
in Fig. 5(b). As for the +A sites, these could have been 
avoided by the trapped particles if the carrier density 
was low. However at n = 0.3, some +A sites also have 
distortions and an associated overdensity, 5n r , but the 
magnitude is smaller than that for —A sites. 
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This happens because locating the localised particles 
always on —A sites would sometimes require occupancy 
of adjacent sites. That would inhibit the (virtual) polaron 
hopping process, and lose kinetic energy. The compro- 
mise is to lose some 'potential energy' and put carriers 
on non adjacent sites which may have unfavourable e». 
This is the origin of the visible 'checkerboard' pattern in 
Fig. 6, where particles generally prefer to remain on the 
diagonal with respect to each other. At stronger disor- 
der, pinning would dominate over kinetic energy consid- 
erations, and +A sites would be depopulated. 

The P(r)) distribution is of course not 5 functions at 
rf^ or r/FL and only has a broad three peak character. 

Based on the data discussed here, as well as the overall 
properties presented earlier 17,20 we suggest the follow- 
ing phenomenology for the T = 0, spin polarised case: 
(a) The weak binary disorder tends to create a landscape 
where electronic eigenfunctions, and the resulting den- 
sity, n r , are weakly inhomogcneous. The strong EP cou- 
pling amplifies this inhomogeneity, by generating lattice 
distortions, leading to a state with a fraction of electronic 
states strongly localised. Figure 7 shows P{x) for several 
A at T — 0.01, and also the inferred density of polaronic 
states n po i at n = 0.3. (a) Polaron formation transfers 
weight in N(u) from around tF to lower energies, creat- 
ing a pseudogap, and reduces the kinetic energy (or effec- 
tive carrier count), as visible in the low frequency optical 
spectral weight 17 n c g(CD,T) = J Q a(u>,T)du>, by a factor 
of 10 (at CD = 1) with respect to the A = probcm. 
(c) There are, however, delocalised but strongly scat- 
tered states which still survive (at moderate A): these 
give rise to finite conductivity at T = 0, and a 'metallic' 
albeit non Drude response in <j(oj). (d) In between the 
strongly localised 'polaronic' states, and the delocalised 
states near the Fermi level, there are possible 'Ander- 
son localised' states as well, by which we mean states 
with large localisation length, £; oc . For 'polaronic' states 
£loc ~ 1- 

A quantitative analysis in terms of these three kinds 
of states is possible, by examining the electronic eigen- 
functions (or inverse participation ratio) in the annealed 
background. We leave such discussion for the future. 

Effects at finite temperature: We have tried to argue 
how the interplay of disorder and EP coupling leads to 
a strongly inhomogeneous state at T = 0. For ease of 
argument let us call this a 'two fluid' scenario, involv- 
ing strongly localised states below the Fermi level, and 
strongly scattered extended states near tp. How does this 
picture evolve at finite T, when the magnetic degrees of 
freedom also come into play? 

The resistivity itself increases quickly with increasing 
temperature, Fig. 1(b), unlike what is observed in the 
non magnetic case 20 where the spin disorder does not 
feed back. We argue that the localised states are now 
subject to additional 'hopping disorder', as evident in 
Fig. 1(a) and Fig. 5(c), and remain ineffective in the 
conduction process. In fact, they are localised marginally 
more due to the additional randomness as the trend in 
the DOS reveals. 
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The extended states, which were earlier scattered by 
the potential fluctuations, r]i, are now also scattered by 
the bond disorder in tij and experience reduced mobil- 
ity. This situation is similar to the previously studied 
interplay of structural disorder and paramagnetic scat- 
tering in the double exchange model, 19,21 where it was 
found that turning on spin disorder indeed increases the 
resistivity of the structurally disordered system. That 
problem, however, dealt with 'quenched' structural and 
magnetic disorder, while both of these are annealed vari- 
ables in the present problem. Furthermore, as the mag- 
netic disorder increases some of extended states (or, more 
accurately, the fraction of extended states in a typical 
equilibrium configuration) reduces. Depending on A and 
A, with increasing T the combined disorder in tij and 
Ci — \xi (see eq.(3)) can drive the 'mobility edge', fi c (T), 
above the Fermi level. This indeed happens at the pa- 
rameter point that we considered in this paper, while 
the full range of possibilities, including a simple 'metal 
to metal' crossover, is detailed in the previous paper. 17 

The high temperature 'insulating' phase that shows 
up for T > Tmit ~ T c is different from a standard 
'Anderson insulator' for two reasons: (i) the localisation 
arises from a complex annealed disorder background, and 
not from an uncorrelated disorder distribution, hence 
the pseudogap features, and (ii) since the thermal in- 
crease in disorder is itself driven by T, the mobility gap 
/Li c (T) — ep is always comparable to T, so there is no 
simple activated behaviour in the high temperature re- 
sistivity. This is roughly consistent with what is observed 
in Lao.7Cao.3Mn03, while the MIT is much stronger, and 
the insulating phase more resistive, in the PrCa systems. 

5. Discussion 

We have focused till now on our own results on the 
MIT obtained in the adiabatic d-H-DE model. It may 
be useful to place these results in relation to previous 
work exploring 'bicriticality' in manganitc models, and 
also comment on cooperative lattice effects and quantum 
fluctuations in spins and phonons. 

(i) The connection with bicriticality: The issue of phase 
competction, the existence of first order phase bound- 
aries, and the effect of disorder in such a regime has been 
explored in simple models. 8,9 The issue was brought to 
focus by the early work of Dagotto and coworkers, 8 who 
suggested that the presence of disorder near a first or- 
der phase boundary between a ferromagnetic metal (FM) 
and an antiferromagnctic (possibly charge ordered) in- 
sulator (AFI) could lead to a pattern of coexisting FM 
and AFI clusters and percolative transport. The detailed 
transport properties in such a system, probed within a 
microscopic theory, were clarified by us. 22 More recently 
a model involving Double Exchange, Holstein EP cou- 
pling, and disorder, has been studied in a two dimen- 
sional system, at half-filling, and discovered the "metalli- 
sation" of an intermediate coupling charge ordered (CO) 
phase by weak disorder. Our own results, far from a CO 
state, suggest how a metal close to a polaronic instability 
is affected by the presence of weak disorder. 

Weak disorder has contrasting roles in these different 
situations, (a) For a generic first order transition, bc- 
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tween a FM and an AFI, say 8 ' 22 disorder acts to con- 
vert "macroscopic" phase separation to meso/nanoscale 
phase coexistence depending on the strength of the ran- 
dom potential. This harks back to the classic Imry-Ma 
scenario 23 and a fascinating variety of percolative effects 
can arise. The key physical effect here is disorder in- 
duced cluster formation near a first order transition and 
does not involve phonons in any essential way. In fact if 
phonons were to be included there are additional effects of 
disorder (beyond Imry-Ma) which need to be considered. 
The second case, 9 (a) involves the effect observed in the 
n = 0.5 CO system. Here disorder generates a pinning 
potential, which acts as a random field, and destroys the 
positional correlations of the CO state. Since the inter- 
mediate coupling CO state depends on the periodicity of 
the CO to generate insulating behaviour, destruction of 
CO promotes metallisation, 9 ' 24 which we have observed 
in our 3d study 20 as well. This is disorder induced met- 
allisation of an intermediate coupling CO state. If the 
EP coupling were large then disorder could still destroy 
the CO phase but the charge ordered insulator would be 
converted to a charge disordered polaronic insulator. Our 
situation, (c) corresponds to a clean metal close to a po- 
laronic instability being converted to a highly resistive 
state by the effect of disorder. Here disorder induces a 
polaronic instability. It should be clear that the physical 
effect in (a) have no essential relation to those in (a) and 
(c). While it is true that both (a) and (c) arise in the same 
model, they occur on quite distinct reference states, and 
the key weak disorder effects are physically very differ- 
ent, despite belonging to the same 'global' A — A — n — T 
phase diagram. 

(ii) Order of the MIT: we think the 'second order' char- 
acter of our MIT is due to the effect of noncooperative 
phonons. If the lattice degrees of freedom were directly 
coupled, as in the real material, then distortions at one 
site can have a cascading effect on the other sites, gen- 
erating an abrupt transition, as indeed observed in one 
MC study. 6 

(iii) Quantum fluctuations in spins and phonons: Our 
T = state has no quantum fluctuations (in spins or 
phonons). While the quantum character of phonons af- 
fects the low T resistivity in clean metals, for the dis- 
ordered strong coupling systems that we consider (with 
large lattice distortions at T = 0) quantum fluctuations 
may not have a qualitative effect. It should be possible 
to make a perturbative estimate of the effects at small 
phonon frequency (u p h/t <C 1), we have not done that 
yet. Similarly, the system is a fully polarised magnet at 
T = 0. Since the effective couplings in the problem are 
ferromagnetic we do not anticipate a significant renor- 
malisation in magnetic properties due to quantum spin 
fluctuations. Finally, by the time T ~ T c , where the MIT 
occurs, the thermal fluctuations in spins and phonons are 
far more important than quantum effects. 

6. Conclusion 

The role of disorder is crucial in the manganites, quite 
independent of its effect at magnetic bicriticality. This is 



due to the inherently strong EP coupling in the mangan- 
ites, which places them close to a polaronic instability. 
Disorder can partially trigger such an instability. The in- 
terplay of disorder and EP coupling controls the resistiv- 
ity and the ferromagnetic T c . Even at fixed density and 
ta (or A) just increasing disorder can completely change 
the transport response. 

While the present investigation has focused only on 
the thermally driven metal-insulator transition, the com- 
bined effect of electron-phonon coupling, disorder, and an 
antiferromagnetic coupling (to compete with double ex- 
change) can generate responses that interpolate between 
the 'canonical DE magnets', Lai_ x Sr x Mn03, to ther- 
mally driven transitions, Lai_ x Ca x Mn03, all the way 
to bicriticality and possible nanoscale coexistence, as in 
Lai-a^PryCaajMnOs. A systematic use of our method 
should clarify the origin, and provide a detailed mi- 
croscopic description, of the wide variety of transport 
regimes observed in the manganites. 

We acknowledge use of the Beowulf cluster at HRI. 
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